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Propagation  of  electromagnetic  waves  through  a  two-dimensional  triangular  lattice  has  been  studied 
for  different  values  of  refractive  index  contrast  between  the  constituent  dielectrics,  and  for  angles  of 
incidence  both  in  and  out  of  the  plane  of  periodicity.  Transmission  results  have  been  obtained  both 
experimentally  and  with  the  transfer  matrix  technique,  and  good  agreement  has  been  found  between 
the  two.  Comparison  with  band  structure  calculations  has  also  been  made.  ©  2001  American 
Institute  of  Physics.  [DOI:  10.1063/1.1330250] 


I.  INTRODUCTION 

Photonic  band  gap  materials  are  artificial  periodic  di¬ 
electric  structures  that  can  control  the  properties  of  light  in  a 
manner  analogous  to  the  way  semiconductor  crystals  control 
the  properties  of  electrons.  Such  dielectric  structures,  some¬ 
times  referred  to  as  “photonic  crystals,”  can  forbid  the 
propagation  of  electromagnetic  (EM)  waves  in  certain  fre¬ 
quency  ranges  called  “photonic  band  gaps,”  which  are  said 
to  be  ‘  ‘complete’  ’  when  the  prohibition  applies  to  all  propa¬ 
gation  directions  and  polarizations.  The  specific  properties  of 
the  band  gaps  (for  example,  their  positions  and  widths)  can 
be  engineered  by  altering  the  geometry  or  lattice  spacing  of 
the  dielectric  structure,  by  adjusting  the  index  contrast  be¬ 
tween  the  constituent  dielectrics,  or  by  introducing  impuri¬ 
ties  or  defects.1'2  The  resulting  versatility  and  unique  prop¬ 
erties  of  photonic  crystals  makes  them  promising  for 
applications1'3  such  as  waveguides,  mirrors,  cavities,  antenna 
substrates,  etc. 

In  previous  work,  it  was  demonstrated  that  complete 
band  gaps  can  exist  in  both  two-  and  three-dimensional  pho¬ 
tonic  crystals  in  specific  cases.4'5  The  two-dimensional  (2D) 
photonic  crystals,  in  particular,  have  attracted  experimental 
interest  because  of  their  advantage  of  simpler  fabrication, 
especially  in  the  case  of  the  small  lattice  spacings  required  to 
produce  band  gaps  at  optical  frequencies.6'7  A  standard  fab¬ 
rication  method  for  air-dielectric  2D  photonic  crystals  is  the 
etching  of  holes  in  a  slab  of  material,  either  by  patterning  a 
dielectric  film  on  a  substrate13  or  by  fabricating  a  channel- 
array  glass.6,9 

It  is  well-known  that  photonic  crystals  with  relatively 
low  refractive  index  contrast  do  not  usually  possess  complete 
photonic  band  gaps.  However,  an  advantage  of  these  materi¬ 
als  is  that  the  directional  stop  bands  occur  at  higher  frequen¬ 
cies  and  are  narrower  than  those  of  corresponding  photonic 
crystals  with  higher  index  contrast.4'610  Additionally,  the  op¬ 


tical  properties  of  low-index-contrast  photonic  crystals  can 
be  useful  for  certain  technological  applications,  such  as  op¬ 
tical  filters2  and  optical  power  limiters.11 

In  a  majority  of  past  studies,  the  optical  properties  of  2D 
photonic  crystals  were  investigated  mainly  for  EM  waves 
propagating  within  the  plane  of  periodicity.  Recently,  how¬ 
ever,  a  lot  of  interest  has  focused  on  the  fabrication  of  optical 
fibers  that  rely  entirely  on  the  existence  of  photonic  band 
gaps  rather  than  on  total  internal  reflection.  In  these  optical 
fibers,  light  is  guided  perpendicularly  to  the  2D  periodic 
plane  through  a  defect  in  the  photonic  crystal  lattice. 12-15 
This  novel  technological  application  of  a  photonic  crystal 
material  suggests  the  need  to  investigate  the  existence  and 
properties  of  photonic  band  gaps  for  light  propagating  out  of 
the  periodic  plane  in  2D  photonic  crystals. 

II.  CALCULATION  METHOD  AND  EXPERIMENT 

In  the  present  work,  experimental  measurements  of  the 
transmission  of  EM  waves  through  a  2D  triangular  photonic 
crystal  are  compared  with  corresponding  theoretical  compu¬ 
tations.  The  method  used  for  the  calculations  is  the  transfer 
matrix  technique  developed  by  J.  B.  Pendry.1617  In  this 
method,  a  grid  is  used  to  divide  the  space  occupied  by  the 
photonic  crystal  into  small  cells.  Maxwell’s  equations  are 
discretized  and  solved  on  this  grid  lattice.  One  is  able  to 
calculate  the  transmission  and  reflection  coefficient  of  a  layer 
consisting  of  only  a  few  grid  points  by  integrating  the  fields 
on  that  layer  along  the  direction  of  propagation.  After  repeat¬ 
ing  this  process  for  adjacent  layers  throughout  the  entire 
sample,  the  total  transmission  coefficient  is  obtained  by  com¬ 
bining  the  coefficients  of  all  the  layers. 

The  experimental  data  were  obtained  by  measuring  the 
transmission  of  channel-glass  photonic  crystals  similar  to 
those  described  previously.-’  ’  -  ’  Measurements  were  per¬ 
formed  along  various  propagation  directions,  both  in  and  out 
of  the  array  plane.  The  real-space  structure  of  these  photonic 
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FIG.  1 .  The  2D  triangular  structure  in  real  space  with  the  corresponding  first 
Brillouin  zone,  showing  normal  incidence  along  (a)  T M  and  (b)  T  K. 


crystals,  along  with  the  corresponding  first  Brillouin  zone 
and  its  respective  high-symmetry  points,  is  depicted  in 
Fig.  1. 

III.  RESULTS  AND  DISCUSSION 

A  series  of  experimental  transmission  spectra  for  several 
in-plane  angles  of  incidence  is  shown  in  Fig.  2.  These  spec¬ 
tra  were  measured  through  a  photonic  crystal  consisting  of 
air  cylinders  in  a  glass  matrix  ( n  —  1.45)  and  having  a  lattice 
constant  of  0.73  /z m.  The  cylinder  diameter  is  approximately 
2/3  of  the  lattice  constant  and  the  cylinder  length  is  approxi¬ 
mately  300  /im.  The  samples  measured  were  approximately 
300  /im  thick  in  the  plane  of  periodicity.  The  spectra  in  Figs. 
2(a)  and  2(b)  were  obtained  with  a  photonic  crystal  whose 
faces  were  oriented  as  shown  in  Fig.  1(a),  at  normal  inci¬ 
dence  and  at  19°  away  from  normal  incidence,  respectively. 
The  spectra  in  Figs.  2(c)  and  2(d)  were  obtained  with  a  pho¬ 
tonic  crystal  whose  faces  were  oriented  as  shown  in  Fig. 
1(b),  at  14°  away  from  normal  incidence  and  at  normal  in¬ 
cidence,  respectively.  Two  polarizations  were  measured  in 
each  case:  electric  field  perpendicular  to  the  cylinder  axis 
(“TE  polarization”),  and  magnetic  field  vector  perpendicu¬ 
lar  to  the  cylinder  axes  (“TM  polarization”).  Similar  experi¬ 
mental  results  for  a  lower  index  contrast  photonic  crystal 
have  been  described  in  Refs.  6  and  10,  and  were  found  to  be 
in  good  agreement  with  band  structure  calculations. 

The  transmission  calculated  via  the  transfer  matrix  tech¬ 
nique  for  the  cases  corresponding  to  Figs.  2(a)-2(d)  is 
shown  in  Figs.  3(a)— 3(d),  respectively.  The  parameters  used 


FIG.  2.  Measured  transmission  for  a  2D  triangular  structure  of  air  cylinders 
in  a  glass  (n~1.45)  matrix.  The  lattice  spacing  is  0.73  /im  and  the  sample 
thickness  is  300  fj. m.  The  results  in  (a)  and  (b)  correspond  to  a  sample 
oriented  as  shown  in  Fig.  1(a),  for  normal  incidence  and  for  incidence  19° 
away  from  the  normal,  respectively.  The  results  in  (c)  and  (d)  correspond  to 
a  sample  oriented  as  shown  in  Fig.  1(b),  for  incidence  at  14°  away  from  the 
normal  and  for  normal  incidence,  respectively.  The  solid  curves  are  for  TM 
polarization  and  the  dotted  curves  are  for  TE  polarization. 


for  the  transfer  matrix  calculations  were  identical  to  those  in 
the  experiments  except  for  cylinder  height  and  the  sample 
thickness.  In  the  calculations,  infinitely  long  cylinders  were 
considered,  and  the  photonic  crystal  consisted  of  only  25  unit 
cells  in  the  direction  normal  to  the  surface.  Thus,  for  a  lattice 
spacing  of  0.73  /im,  the  calculations  effectively  corre¬ 
sponded  to  a  thickness  of  —  31  /im  when  the  photonic  crystal 
was  oriented  as  shown  in  Fig.  1(a),  and  to  —18  /xm  when 
oriented  as  shown  in  Fig.  1(b).  However,  it  should  be  noted 
that  the  sample  thickness  does  not  affect  the  position  of  the 
photonic  band  gap  but  only  the  strength  of  the  attenuation.  A 
related  issue  is  that  of  sample  termination,  i.e.,  the  positions 
of  the  sample  faces  relative  to  the  centers  of  the  nearest  rows 
of  cylinders.  Although  the  sample  termination  could  not  be 
controlled  experimentally  due  to  the  sub-micron  spacing  be¬ 
tween  cylinders,  in  the  calculations  each  sample  face  was 
assumed  to  occur  near  the  center  of  a  row  of  cylinders,  for 
convenience.  Finally,  it  should  also  be  noted  that  an  “al- 
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Wavelength  (nm) 

FIG.  3.  Transmission  results  for  in-plane  propagation  obtained  with  the 
transfer  matrix  technique.  The  curves  in  (a)-(d)  correspond  to  the  results  in 
Figs.  2(a)— 2(d),  respectively.  The  structure  consists  of  25  unit  cells  in  the 
I’M  or  \  K  direction  when  the  sample  is  oriented  as  shown  in  Fig.  1(a)  or 
Fig.  1(b),  respectively.  The  solid  curves  are  for  TM  polarization  and  the 
dotted  curves  are  for  TE  polarization. 

most’  ’  triangular  structure  was  considered  in  the  calculations 
instead  of  a  perfectly  triangular  one,  for  reasons  that  will  be 
explained  in  detail  below. 

Before  proceeding  with  the  comparison  of  the  measured 
and  calculated  transmission  results,  it  is  worthwhile  to  com¬ 
pare  the  present  calculated  results  with  the  photonic  band 
structure  as  obtained  by  the  plane  wave  expansion  method.4'5 
The  same  photonic  crystal  parameters,  listed  above,  are  used 
for  both  calculations.  For  TE  polarization,  the  band  structure 
calculations  predict  a  band  gap  from  —1466  nm  to  —1710 
nm  at  the  M  symmetry  point,  and  from  — 1315  nm  to  —  1490 
nm  at  the  K  symmetry  point.  For  TM  polarization,  the  band 
structure  calculations  predict  a  band  gap  from  —  1567  nm  to 
— 1738  nm  at  the  M  symmetry  point,  and  no  band  gap  at  the 
K  symmetry  point. 

A  small  discrepancy  between  the  results  of  the  transmis¬ 
sion  and  band  structure  calculations  arises  from  the  limita¬ 
tions  of  the  transfer  matrix  method.  For  this  calculational 
method,  a  finite  number  of  square  cells  (of  the  grid  lattice 
used  to  discretize  Maxwell’s  equations)  is  required  to  fit 
within  a  unit  cell  of  the  2D  photonic  crystal.  Thus,  the  trian¬ 
gular  lattice,  which  is  essentially  a  centered  rectangular  lat¬ 
tice  with  a  lattice  ratio  of  \/3,  has  to  be  approximated  by  one 
with  a  lattice  ratio  of  1.7.  This  would  be  responsible  for  a 
discrepancy  up  to  2% .  Also,  the  finite  (5  nm)  step  taken  for 
the  wavelength  when  calculating  the  transmission  imposes 
an  error  on  the  value  of  the  band  gap  edges  obtained  from 
the  transmission  results. 

For  TE  polarization,  the  transmission  calculated  by  the 
transfer  matrix  method  determines  the  lowest-energy  band 
gap  to  be  from  —1440  nm  to  —1682  nm  for  normal  inci¬ 


dence  along  T M,  as  shown  in  Fig.  3(a),  and  from  —  1325  nm 
to  —1480  nm  for  normal  incidence  along  F  K ,  as  shown  in 
Fig.  3(d).  These  results  agree  within  5%  with  those  of  the 
band  structure  calculation.  However,  the  case  of  TM  polar¬ 
ization  is  more  complex.  Although  the  transmission  calcula¬ 
tion  determines  a  band  gap  from  —1530  nm  to  —1700  nm 
for  normal  incidence  along  FM,  as  shown  in  Fig.  3(a), 
which  is  also  within  5%  of  the  result  of  the  band  structure 
calculation,  the  transmission  calculation  also  predicts  a  band 
gap  from  — 1235  nm  to  — 1456  nm  for  incidence  along  T  K , 
as  shown  in  Fig.  3(d),  where  the  band  structure  predicts  no 

gap- 

Experimentally,  for  TE  polarization,  the  measured  band 
gaps  appear  from  — 1335  nm  to  — 1559  nm  for  normal  inci¬ 
dence  along  FM  in  Fig.  2(a),  and  from  —  1 154  nm  to  —  1352 
nm  for  normal  incidence  along  \  K  in  Fig.  2(d).  For  TM 
polarization,  the  measured  band  gaps  appear  from  — 1422 
nm  to  —1584  nm  for  normal  incidence  along  FM  in  Fig. 
2(a),  and  from  —1127  nm  to  —1358  nm  for  normal  inci¬ 
dence  along  FK  in  Fig.  2(d).  This  represents  a  discrepancy 
between  the  measured  and  calculated  transmission,  described 
above,  ranging  from  7%  to  13%.  Similar  agreement  is  also 
found  for  the  cases  of  off-normal  incidence,  i.e.,  comparing 
Figs.  2(b)  and  2(c)  to  Figs.  3(b)  and  3(c),  respectively. 
Known  sources  of  experimental  errors  include  uncertainties 
in  the  angle  of  incidence  (±1  degree),  in  the  lattice  spacing 
(±5%),  and  in  the  glass  index  of  refraction  (±5%),  which 
is  modified  by  the  etching  process.  Since  the  position  of  the 
band  gap  is  proportional  to  the  lattice  spacing,  errors  in  the 
lattice  spacing  are  linearly  proportional  to  errors  in  the  band 
gap  position.  Also,  from  a  transfer  matrix  calculation,  it  was 
found  that  a  3%  change  in  the  glass  index  can  cause  a  cor¬ 
responding  change  of  up  to  3%  in  the  positions  of  the  band 
edges.  These  experimental  uncertainties,  in  combination  with 
the  limitations  of  the  transfer  matrix  method  discussed 
above,  can  explain  the  discrepancies  between  the  transmis¬ 
sion  spectra  measured  experimentally  and  those  calculated 
with  the  transfer  matrix  method. 

Nevertheless,  the  transmission  (obtained  both  experi¬ 
mentally  and  through  transfer  matrix  calculations)  remains 
inconsistent  with  the  band  structure  results  by  showing  a  gap 
at  the  K  symmetry  point  where  the  band  structure  predicts 
none.  This  type  of  apparent  contradiction  between  the  band 
structure  results  and  the  transmission  (measured  and  calcu¬ 
lated)  of  a  photonic  crystal  has  been  observed  in  previous 
experiments.19'20  In  Refs.  19  and  20,  it  is  shown  that  even 
when  a  certain  photonic  crystal  can  host  a  particular  mode, 
that  mode  will  not  be  excited  by  the  incident  field  if  the  field 
patterns  which  correspond  to  that  mode  have  symmetry  in¬ 
consistent  with  that  of  the  incident  field.  This  causes  a  pseu¬ 
dogap  to  appear  in  the  transmission.  Also,  Sakoda2122  has 
performed  calculations  for  a  triangular  lattice  of  air  cylinders 
in  a  glass  matrix  ( n~  1.65).  He  showed  that  the  second  band 
for  the  TM  polarization  (which  is  degenerate  at  the  K-point 
along  the  F  K  direction)  will  not  couple  at  all  to  an  incoming 
field  incident  along  the  F K  direction  when  that  direction  is  at 
normal  incidence  for  the  photonic  crystal.  Taking  this  into 
consideration,  one  expects  the  appearance  of  a  pseudogap  in 
the  photonic  crystals  used  in  Fig.  2,  between  — 1263  nm  and 
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FIG.  4.  Measured  transmission  results  for  a  2D  triangular  structure  of  air 
cylinders  in  a  glass  (n~  1.45)  matrix.  The  propagation  direction  is  at  45° 
out  of  the  array  plane  with  respect  to  TM  in  (a),  and  within  the  array  plane 
along  TM  in  (b).  Two  samples  were  used,  each  having  the  propagation 
direction  at  normal  incidence.  The  lattice  spacing  is  0.56  /zm  and  the  sample 
thickness  is  300  /zm  in  each  case.  The  solid  curves  are  for  TM  polarization 
and  the  dotted  curves  are  for  TE  polarization. 


~  1496  nm  for  TM  polarization  and  normal  incidence  along 
the  TAT  direction.  This  is  in  very  good  agreement  with  the 
transmission  results  (obtained  both  experimentally  and 
through  transfer  matrix  calculations),  and  also  explains  why 
for  this  incident  direction  the  TM  polarization  gap  appears  to 
be  wider  than  the  TE  polarization  gap. 

In  the  calculated  TM  polarization  transmission  spectrum 
of  Fig.  3(c)  there  is  another  striking  feature  which  deserves 
further  discussion.  Besides  a  gap  from  ~  1475  nm  to  ~  1620 
nm  there  is  a  “dip”  in  the  transmission  for  the  wavelength 
region  between  —1130  nm  and  — 1295  nm.  This  feature  cor¬ 
relates  well  with  a  similar  feature  observed  experimentally  in 
Fig.  2(c).  The  transmission  in  the  wavelength  region  of  the 
dip  differs  significantly  from  zero  and  also  shows  interfer¬ 
ence  fringes.  By  contrast,  the  transmission  in  the  gap  region 
for  a  sufficiently  thick  sample  (as  is  the  case  for  the  sample 
used  for  the  calculations)  is  always  nearly  zero.  Additional 
transfer  matrix  calculations  have  shown  that  the  attenuation 
associated  with  the  “dip”  is  independent  of  sample  thick¬ 
ness,  while  the  attenuation  associated  with  either  a  true  gap 


FIG.  5.  Measured  transmission  results  similar  to  Fig.  4,  but  for  a  structure 
consisting  of  water-filled  (w  —  1 .33)  cylinders  in  a  glass  (/r  —  1 .45)  matrix. 
The  propagation  direction  is  at  45°  out  of  the  array  plane  with  respect  to 
TM  in  (a),  and  within  the  array  plane  along  TM  in  (b).  The  additional 
absorption  appearing  near  1450  nm  in  both  (a)  and  (b)  is  due  to  water.  The 
solid  curves  are  for  TM  polarization  and  the  dotted  curves  are  for  TE  po¬ 
larization. 


or  of  a  “pseudogap”  [like  the  one  in  Fig.  3(d)]  does  depend 
on  sample  thickness.  Also,  further  calculations  have  shown 
that  the  fringes  within  the  “dip”  as  well  as  the  magnitude  of 
the  attenuation  in  this  spectral  region  are  both  sensitive  to  the 
sample  termination.  In  addition,  the  band  structure  calcula¬ 
tions  show  that  propagating  modes  do  exist  in  the  wave¬ 
length  region  of  the  “dip.”  These  observations  lead  to  the 
conclusion  that  the  observed  transmission  “dip”  is  likely  to 
result  from  the  “partial”  (i.e.,  incomplete)  coupling  between 
the  incident  wave  and  the  propagating  modes.  (Note  that,  if 
the  coupling  were  complete,  one  would  observe  ordinary 
transmission  fringes  oscillating  up  to  the  maximum  transmis¬ 
sion  value  of  1.) 

The  experimental  transmission  spectra  for  incidence  di¬ 
rections  outside  the  array  plane  for  similar  2D  triangular 
photonic  crystals  with  several  index  contrasts  are  shown  in 
Figs.  4,  5,  and  6.  The  measured  in-plane  transmission  corre¬ 
sponding  to  each  index  contrast  is  also  included  for  refer¬ 
ence.  The  photonic  crystals  used  for  these  measurements 
were  also  based  on  channel  glass  materials,  as  above,  except 
oriented  so  that  the  incident  propagation  direction  of  interest 
was  along  the  surface  normal  for  each  sample.  The  particular 
direction  chosen  for  the  data  in  Figs.  4(a),  5(a)  and  6(a)  is 
45°  away  from  the  in-plane  F M  direction.  The  lattice  spac¬ 
ing  of  these  photonic  crystal  samples  was  0.56  /im,  and  the 
thickness  of  each  sample  was  —300  /rm.  Results  for  the 
highest  index  contrast  achievable,  namely  for  a  photonic 
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FIG.  6.  Measured  transmission  results  similar  to  Fig.  4,  but  for  a  structure 
consisting  of  glass  («~1.67)  cylinders  in  a  glass  (« —  1 .65)  matrix.  The 
propagation  direction  is  at  45°  out  of  the  array  plane  with  respect  to  TM  in 
(a),  and  within  the  array  plane  along  TM  in  (b).  The  solid  curves  are  for  TM 
polarization  and  the  dotted  curves  are  for  TE  polarization. 


crystal  consisting  of  air  cylinders  in  a  glass  (n  —  1.45)  ma¬ 
trix,  are  shown  in  Fig.  4.  An  intermediate  index  contrast  was 
obtained  in  Fig.  5  by  immersing  the  photonic  crystal  from 
Fig.  4  in  water  (n~  1.33).  (Additional  absorptions  due  to  the 
water,  which  overlap  slightly  with  the  attenuations  associ¬ 
ated  with  the  photonic  band  gaps,  appear  in  Fig.  5).  The 
lowest  index  contrast  (Fig.  6)  was  obtained  in  a  photonic 
crystal  consisting  of  glass  («  —  1 .67)  cylinders  in  a  glass  ( n 
■ —  1 .65)  matrix.  (Note  that  all  glass  refractive  indices  are 
wavelength-dependent;  the  indices  quoted,  which  were  used 
in  the  calculations,  are  accurate  at  —800  nm).  Also  note  that 
a  more  detailed  study  of  the  out-of-plane  transmission  for  the 
lowest-index-contrast  case  has  been  published  in  Ref.  18. 

The  calculated  transmission  spectra  corresponding  to 
Figs.  4  and  6,  obtained  via  the  transfer  matrix  technique  for 
incidence  in  the  array  plane  and  at  a  45°  angle  outside  the 
array  plane  with  respect  to  the  in-plane  TM  direction,  are 
shown  in  Figs.  7  and  8.  Although  some  fine  structure  is 
associated  with  the  attenuations  appearing  in  the  transmis¬ 
sion  calculations  for  the  lowest  index  contrast  photonic  crys¬ 
tal,  at  —810  nm  for  in-plane  incidence  and  at  —725  nm  for 
out-of-plane  incidence,  only  one  (broader)  attenuation  is 
measured  experimentally.  This  is  at  least  partially  a  result  of 
the  slight  disorder  inevitably  present  in  the  actual  photonic 
crystals,  which  results  in  inhomogeneous  broadening  of 
these  narrow  attenuation  regions.  Therefore,  for  comparison 
with  the  experimental  results,  only  the  center  values  of  the 
calculated  attenuation  features  were  considered.  Very  good 


Wavelength  (nm) 


FIG.  7.  Transmission  results  for  out-of-plane  propagation  obtained  with  the 
transfer  matrix  technique,  for  a  triangular  photonic  crystal  consisting  of  air 
cylinders  in  a  glass  matrix,  corresponding  to  Fig.  4.  The  curves  in  (a)  and  (b) 
correspond  to  the  results  in  Figs.  4(a)  and  4(b),  respectively.  Note  that  the 
boundary  conditions  for  (a)  differ  from  those  of  Fig.  4(a).  The  sample  sur¬ 
face  is  at  45°  with  the  incident  wave  vector  while  in  Fig.  4(a)  the  surface  is 
perpendicular  to  the  wave  vector.  The  solid  curves  are  for  TM  polarization 
and  the  dotted  curves  are  for  TE  polarization. 

agreement  is  found  between  the  measured  transmission  and 
the  transfer  matrix  results.  The  discrepancies  range  from 
—  1  %  to  6%  for  in-plane  incidence  and  from  —6%  to  9%  for 
out-of-plane  incidence.  The  larger  disagreement  in  the  out- 
of-plane  case  is  expected:  The  cylinders  in  the  experimental 
sample  had  finite  length  and,  for  the  measurements  presented 
in  Figs.  4(a),  5(a),  and  6(a),  they  were  cut  in  such  a  way  that 
incidence  is  normal  to  the  sample  surface  in  each  case,  as 
mentioned.  By  contrast,  the  numerical  method  considered 
infinitely  long  cylinders  that  are  parallel  to  the  sample  sur¬ 
face,  which  is  oriented  at  45°  with  respect  to  the  incident 
wave  vector.  Since  the  boundary  conditions  imposed  on  the 
sample  surface  require  that  the  parallel  component  of  the 
incident  wave  vector  be  conserved  in  the  crystal,23  the  result¬ 
ing  transmission  can  be  influenced  by  the  orientation  of  the 
sample  surface  relative  to  the  incident  wave  vector.  In  spite 
of  these  differences,  however,  good  agreement  was  achieved 
between  the  transfer  matrix  calculations  and  the  experimen¬ 
tal  results. 

An  important  observation  is  that  the  photonic  band  gaps 
(or  the  corresponding  attenuations)  for  all  index  contrast  val¬ 
ues  persist  for  incidence  out  of  the  array  plane  but  shift  to¬ 
ward  shorter  wavelengths.  Flowever,  the  depths  of  all  of  the 
attenuations  appear  to  decrease  in  the  measured  spectra. 
Such  a  decrease  was  not  found  in  the  transfer  matrix  results, 
except  for  the  case  of  the  lowest  index  contrast  (correspond¬ 
ing  to  glass  cylinders  in  a  glass  matrix),  and  then  only  for  the 
lowest-frequency  attenuation  in  the  TM  polarization.  (The 
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FIG.  8.  Transmission  results  for  out-of-plane  propagation  obtained  with  the 
transfer  matrix  technique,  for  a  triangular  photonic  crystal  consisting  of 
glass  cylinders  in  a  glass  matrix,  corresponding  to  Fig.  6.  The  curves  in  (a) 
and  (b)  correspond  to  the  results  in  Figs.  6(a)  and  6(b),  respectively.  Note 
that  the  boundary  conditions  for  (a)  differ  from  those  of  Fig.  6(a).  The 
sample  surface  is  at  45°  with  the  incident  wave  vector  while  in  Fig  6(a)  the 
surface  is  perpendicular  to  the  wave  vector.  The  solid  curves  are  for  TM 
polarization  and  the  dotted  curves  are  for  TE  polarization.  In  (b)  the  dotted 
curve  is  covered  by  the  solid  curve. 

angular  dependence  of  the  strength  of  this  attenuation  was 
discussed  in  detail  in  Ref.  18).  Another  important  point  is 
that  calculations  for  the  highest-index-contrast  case  (corre¬ 
sponding  to  air  cylinders  in  a  glass  matrix),  performed  for 
various  lattice  filling  ratios,  show  that  the  widest  gaps  occur 
when  the  filling  ratio  is  between  40%  and  50%  (defined  as 
the  sum  of  the  cross-sectional  areas  of  the  cylinders  over 
total  area).  The  experiments  presented  here  were  done  for  a 
cylinder  diameter  approximately  2/3  of  the  lattice  spacing, 
which  corresponds  to  a  40%  filling  ratio.  Thus,  for  the 
present  case  of  the  triangular  lattice  of  air  cylinders  in  a  glass 
matrix,  the  optimal  band  gap  is  obtained. 

IV.  CONCLUSIONS 

The  transmission  of  2D  triangular  lattices,  with  index 
contrast  corresponding  to  air  cylinders  in  a  glass  matrix  and 
below,  was  studied  both  experimentally  and  through  transfer 
matrix  calculations.  The  experimental  results  show  generally 
good  agreement  with  the  calculations,  with  discrepancies 
ranging  from  1%  to  13%.  A  narrow  band  gap  region  inde¬ 
pendent  of  incidence  direction  was  found  only  for  the  pho¬ 
tonic  crystals  consisting  of  air  cylinders  in  a  glass  matrix, 
and  only  for  in-plane  propagation  and  TE  polarization.  None 
of  the  dielectric  structures  considered  here  possesses  a  com¬ 
plete  2D  photonic  band  gap  (independent  of  polarization,  for 
in-plane  propagation).  For  TM  polarization,  the  appearance 
of  a  ‘  ‘pseudogap’  ’  was  observed  both  in  the  measured  spec¬ 


tra  and  in  the  transfer  matrix  calculations.  The  positions  of 
the  band  gaps  corresponding  to  the  lowest  index  contrast 
structures,  consisting  of  glass  cylinders  in  a  glass  matrix, 
depend  on  the  incident  direction  and  are  relatively  narrow, 
making  them  useful  for  applications  such  as  optical  filters.2 
The  photonic  crystals  studied  in  this  paper  retain  their  band 
gaps  even  for  out-of-plane  propagation,  but  the  gaps  are 
shifted  toward  shorter  wavelengths.  The  absence  of  overlap 
between  the  out-of-plane  and  the  in-plane  band  gaps,  even  in 
the  case  of  the  photonic  crystals  consisting  of  air  cylinders  in 
a  glass  matrix,  suggests  that  this  particular  structure  is  inap¬ 
propriate  for  use  as  a  photonic  band  gap  optical  fiber.  How¬ 
ever,  previous  transfer  matrix  calculations24  for  a  triangular 
photonic  crystal  with  a  dielectric  contrast  of  16  have  shown 
the  existence  of  a  complete  band  gap  for  incidence  up  to  85° 
out  of  the  array  plane.  Furthermore,  understanding  the  propa¬ 
gation  of  EM  waves  at  arbitrary  directions  through  a  photo¬ 
nic  crystal  should  be  a  useful  step  toward  incorporating 
finite-size  photonic-crystal  elements  into  the  next  generation 
of  opto-electronic  circuits. 
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